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lateral  and  longitudinal  swashplate  angle,  respectively,  deg 
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acceleration  along  body  X  and  Y  axis,  respectively,  as  measured 
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SUMMARY 


■  nJ 

This  report  documents  a  mathematical  model  of  an  advanced  helicopter;  the 
model  is  suitable  for  use  in  control/display  research  involving  piloted  simu¬ 
lation.  The  general  design  approach  for  the  six-degree-of-freedom  equations 
of  motion  is  to  use  the  full  set  of  nonlinear  gravitational  and  inertial  terms 
of  the  equations  and  to  express  the  aerodynamic  forces  and  moments  as  the 
reference  values  and  first-order  terms  of  a  Taylor  series  expansion  about  a 
reference  trajectory  defined  as  a  function  of  longitudinal  airspeed.  Provi¬ 
sions  for  several  different  specific  and  generic  flight  control  systems  are 
included  in  the  model.  The  logic  required  to  drive  various  flight  control  and 
weapon  delivery  symbols  on  a  pilot's  electronic  display  is  also  provided. 
Finally,  the  model  includes  a  simplified  representation  of  low-altitude  wind 
and  turbulence  effects.  This  model  has  been  used  in  a  piloted  simulator 
investigation,  recently  conducted  at  Ames  Research  Center,  of  the  effects  of 
control-system  and  display  variations  for  an  attack  helicopter  mission. 
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1 .  INTRODUCTION 


To  provide  a  real-time  simulation  for  research  programs  that  involve 
wide  variations  in  helicopter  flight  control  system  characteristics  and  elec¬ 
tronic  display  format  and  logic,  it  is  necessary  to  produce  a  model  of  the 
helicopter  that  requires  a  minimum  of  digital  computer  capacity  yet  is  suffi¬ 
ciently  realistic  for  the  task  being  simulated.  One  candidate  technique  has 
been  applied  to  two  fixed-base  simulations  for  control/display  research  involv¬ 
ing  VTOL  aircraft  in  terminal-area  operations  (refs.  1,  2);  this  approach 
includes  the  full  nonlinear  set  of  kinematic  terms  in  six-degree-of-freedom 
equations  of  motion,  but  uses  linearized  aerodynamic  terms  assuming  small 
perturbations  about  reference  conditions  that  vary  as  functions  of  a  single 
independent  variable.  As  described  herein,  this  particular  technique  was 
applied  to  derive  a  mathematical  model  of  an  advanced  helicopter  valid  for 
small  perturbations  from  level  flight  conditions  for  the  airspeed  range  of 
-40  to  +160  knots.  Sufficient  flexibility  is  retained  in  the  model  to  facili¬ 
tate  the  conversion  to  different  basic  helicopter  aerodynamic  characteristics. 


J 

i 

ii 


Bgg 


-“1 


The  derivation  of  the  expressions  for  the  aerodynamic  forces  and  moments  is 
presented  in  section  2  of  this  report. 

A  separate  segment  of  the  model  is  used  to  generate  the  control  surface 
positions  as  functions  of  cockpit  control  positions  and  elements  of  the  aircraft 
state  vector.  A  total  of  six  different  flight  control  systems  is  provided. 

Their  characteristics  are  described  in  section  3. 

Another  portion  of  the  model  contains  the  logic  that  converts  the  elements 
of  the  aircraft  state  vector  into  the  parameters  required  to  drive  the  symbols 
on  the  pilot's  electronic  display.  These  parameters  —  which  convey  aircraft 
orientation,  situation,  command,  and  fire  control  information  to  the  pilot  — 
are  defined  in  section  4. 

Section  5  describes  the  model  of  the  low-altitude  wind  and  turbulence 
implemented  to  provide  external  disturbances  to  the  systems  under  investiga¬ 
tion.  The  model  includes  steady  wind,  wind  shear  in  magnitude  and  direction, 
and  turbulence  components  whose  characteristics  vary  with  mean  wind  speed  and 
altitude.  Conclusions  are  presented  in  section  6. 

These  elements  of  the  mathematical  model  of  an  advanced  helicopter  have 
been  used  in  a  piloted  moving-base  simulator  study  of  the  effects  on  flying 
qualities  of  variations  in  both  flight  control  system  and  display  characteris¬ 
tics  for  an  attack  helicopter  mission  (ref.  3). 


2 .  AERODYNAMICS 


General 

The  total  aerodynamic  forces  and  moments  required  for  the  six-degree-of- 
freedom  equations  of  motion  are  generated  as  the  summation  of  reference  and 
first-order  terms  of  a  Taylor  series  expansion  about  a  reference  trajectory 
defined  as  a  function  of  longitudinal  airspeed  (ua).  Function  generation  sys¬ 
tem  subroutines  (ref.  4)  are  utilized  to  produce  the  values  for  the  following 
parameters  as  functions  of  the  single  variable  ua: 

1.  Reference  values  for  total  forces  and  moments  —  XR,  YR,  ZR,  and  MR 

2.  Reference  values  for  aircraft  motion  and  control  variables  — 

wr’  AV  BV  ev  and  H 

3.  Values  for  the  aircraft  stability  and  control  parameters  —  e.g., 

and  Zp 
°o 

The  reference  values  for  the  total  forces  and  moments  are  specified  at 
20-knot  intervals  of  the  independent  variable  for  -40  knots  <  ua  <  +160  knots. 
Each  of  the  remaining  dependent  variables  is  specified  at  40-knot  intervals  of 
the  independent  variable.  Linear  interpolation  is  used  to  determine  the  value 
of  each  parameter  between  these  breakpoints. 
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Approach 

As  an  example  of  the  technique  used  to  calculate  the  total  aerodynamic 
forces  and  moments ,  the  X-force  equation  is  expressed  as: 


—  =  XR(ua)  +  Xw(ua)Awa  +  Xq(ua)q  +  X^  (ua)ABls  +  X0  (ua)A0o 


where 


Awa  -  wa  -  w 


(u_) 


ABls  =  Bls  '  BlsR(ua) 

A0o  =  eo  “  9oD(ua> 

R 


Preliminary  computer- generated  values  of  the  stability  and  control 
derivatives  and  reference  values  of  the  aircraft  motion  and  control  variables 
have  been  provided  by  Hughes  Helicopters  as  a  function  of  flight  condition; 
values  derived  from  the  flight-test  data  presented  in  reference  5  supplement 
these  computer-generated  data  where  possible. 


The  reference  value  of  the  X-force  (XR)  is  calculated  by  taking  the  par¬ 
tial  derivative  of  Xa  with  respect  to  ua: 


m  3ua  -  ^  dua  dua  ABi  dua 


d0. 


-  X 


R 


eo  dua 


Transposing, 


dX  dw 

=  X,,  +  X,,  -3-^  +  X„ 
dua  u  ^  dua  Bxs 


dBi 


d0. 


+  X 


R 


0O  dua 


Therefore,  using  the  dummy  variable  a: 


XR(ua)  =  XR(ua  =  0)  + 


r 


do 


do 


where  XR(ua  =  0)  =  g  sin  0R  and  0R  =  hover  trim  pitch  attitude. 

For  a  given  airspeed,  the  accuracy  of  the  linear  approximations  of  the 
total  aerodynamic  forces  and  moments  is  dependent  in  part  on  the  magnitude  of 
the  perturbations  from  the  reference  trajectory.  No  analysis  was  conducted  to 
investigate  the  errors  that  result  from  neglecting  the  nonlinear  effects  caused 
by  large  perturbations.  Therefore,  to  improve  confidence  in  the  validity  of 
experimental  resuKs,  the  task  to  be  performed  should  in  general  involve  only 
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"small"  perturbations  from  the  selected  reference  trajectory;  this  is  particu¬ 
larly  true  when  the  model  is  used  for  piloted  simulation. 


Ground  Effect 


To  account  for  the  effect  of  the  terrain  on  the  aircraft  dynamics,  the 
stability  parameter  is  included  in  the  model.  The  derivation  of  the 

value  of  this  parameter  proceeds  as  follows.  From  momentum  theory: 


2C„ 


ao 


For  y  =  0: 


For  the  advanced  helicopter. 


From  reference  6,  the  value  of  X 


CT  =  0.0066 
a  -  5.75 
a  «  0.092 

in  ground  effect  is  expressed  as: 
V  -  30^ 


X  -  »o[l  +  <0  -  1)  (So22)  ] 


For  V  =  0, 


where 


X  =  GX0 


-0.0575 


G  =  0.25  + 


eg 

48 


Therefore , 


and 


0o  =  0.09653  +  0.001797  hcg,  rad 

39o 

-TT-  -  0.001797,  rad/ft 
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Therefore,  for  V  =  0: 


Z,  =  -Za  -  (4.562)  (0.103)  =  0.47  ft/sec2/ft 

n  u0  dn 

According  to  reference  6,  this  effect  diminishes  to  zero  as  airspeed 
approaches  30  ft/sec  and  is  only  apparent  for  altitudes  less  than  1  rotor 
diameter.  In  the  model,  this  effect  is  only  present  for  altitudes  less  than 
50  ft  AGL.  (Because  no  direct  measurement  of  height  above  the  simulated  ter¬ 
rain  is  available,  altitude  AGL  is  estimated  by  subtracting  a  constant  value 
representing  mean  terrain  altitude  from  the  actual  altitude  above  the  terrain 
board  reference.)  The  parameter  decreases  linearly  with  velocity,  reach¬ 

ing  a  value  of  zero  for  ua  =  -40  and  +40  knots. 


Coupling  Derivatives 

The  longitudinal  and  lateral-directional  aerodynamics  of  the  basic  model 
are  uncoupled  with  the  exception  of  yawing  moment  due  to  tail  rotor  collective 
pitch  inputs.  An  option  which  adds  perturbations  to  the  basic  aerodynamic 
forces  and  moments  to  account  for  coupling  effects  is  available.  The  follow¬ 
ing  coupling  effects  are  included:  (1)  longitudinal  equations,  vfl,  p,  r, 

Ais,  0TR;  and  (2)  lateral-directional  equations,  wa,  q,  0Q,  Bis. 


Summary  of  Equations 

Perturbation  variables- 

DWB  =  WB  -  WBR 
DTHET0  =  THET0  -  THET0R 
DAIS  =  AIS  -  AlSR 
DB1S  =  BIS  -  B1SR 
DTHETTR  =  THETTR  -  THETTRR 

where  WBR,  THET0R,  AlSR,  B1SR,  and  THETTRR  are  all  generated  by  function 
generator  system  subroutines  as  functions  of  UB. 

X-force  equation- 

FAX  -  XMASS*{XQ*QB  +  XW*DWB  +  XB1S*DB1S  +  XTH0*DTHET0  +  XREF} 


where  XQ,  XW,  XB1S,  XTH0,  and  XREF  are  all  generated  as  functions  of  UB. 


Y-force  eguation- 

FAY  =  XMASS*{YP*PB  +  YR*RB  +  YV*VB  +  YA1S*DA1S  +  YTHTR*DTHETTR  +  YREF ; 
where  YP,  YR,  YV,  YA1S,  YTHTR,  and  YREF  are  all  generated  as  functions  of  UB. 
Z-force  equation- 

FAZ  =  XMASS*{ZQ*QB  +  ZW*DWB  +  ZB1S*DB1S  +  ZTH0*DTHET0  +  ZH*DH  +  ZREF } 


where  ZQ,  ZW,  ZB1S,  ZTH0,  ZH,  and  ZREF  are  all  generated  as  functions  of  UB 
and 


DH 


HAGL  -50  for  HAGL  <  50  ft 
0  for  HAGL  >  50  ft 


where  HAGL  =  HCG  -  HTER 


L-moment  equation- 

TAL  =  XIXX* { ULP*PB  +  ULR*RB  +  ULV*VB  +  ULA1S*DA1S  +  ULTTR*DTHETTR } 


where  ULP,  ULR,  ULV,  ULA1S,  and  ULTTR  are  all  generated  as  functions  of  UB. 
M-moment  equation- 

TAM  =  XIYY* { UMQ*QB  +  UMW*DWB  +  UMB1S*DB1S  +  UMTH0*DTHET0  +  UMREF } 
where  UMQ,  UMW,  UMB1S,  UMTH0,  and  UMREF  are  all  generated  as  functions  of  UB. 
N-moment  equation- 

TAN  =  XIZZ*{UNP*PB  +  UNR*RB  +  UNV*VB  +  UNTH0*DTHET0 
+  UNTTR*DTHETTR  +  UNA1S*DA1S} 

where  UNP,  UNR,  UNV,  UNTH0,  UNTTR,  and  UNA1S  are  all  generated  as  functions  of 
UB. 


The  values  of  the  reference  forces  and  moments,  stability  and  control 
parameters,  and  reference  aircraft  motion  and  control  variables  are  presented 
in  the  appendix  as  functions  of  ua  at  the  designated  breakpoints. 


The  optional  perturbations  to  the  basic  expressions  for  total  aerodynamic 
forces  and  moments  to  account  for  coupling  effects  are: 

DELFAX  =  XMASS*{UKPBX*PB  +  UKRBX*RB  +  UKVBX*VB  +  UKDA1SX*DA1S 

+  UKDTHTRX*DTHETTR } 
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DELFAY  =  XMASS*{UKQBY*QB  +  UKDWBY*DWB  +  UKDB1SY*DB1S  +  UKDTH0Y  *DTHET0 } 

DELFAZ  =  XMASS*{UKPBZ*PB  +  UKRBZ*RB  +  UKVBZ*VB  +  UKDA1SZ*DA1S 
+  UKDTHTRZ*DTHETTR } 

DELTAL  =  XIXX*{UKQBL*QB  +  UKDWBL*DWB  +  UKDB1SL*DB1S  +  UKDTH0L*DTHET0 } 

DELTAM  =  XIYY* { UKPBM*PB  +  UKRBM*RB  +  UKVBM*VB  +  UKDA1SM*DA1S 
+  UKDTHTRM*DTHETTR} 

DELTAN  =  XIZZ*{UKQBN*QB  +  UKDWBN*DWB  +  UKDB1SN*DB1S  +  UKDTH0N*DTHET0 } 

The  values  for  the  coupling  derivatives  are  presented  in  the  appendix. 

3.  FLIGHT  CONTROL  SYSTEMS 

General 

Six  control  system  configurations  are  provided: 

1.  Mechanical  —  pitch,  roll,  yaw  SCAS  OFF 

2.  SCAS  ON  —  pitch,  roll,  yaw;  individually  or  collectively 

3.  Attitude  Hold  —  pitch,  roll,  yaw 

4.  Hover  augmentation  —  pitch  and  roll 

5.  Heading  hold  —  yaw 

6.  Vertical  augmentation  —  collective 

Configurations  (1)  through  (3)  are  specific  control  system  configurations 
based  in  part  on  information  derived  from  references  5  and  7  and  on  unpublished 
stability  and  control  data.  Configurations  (4)  through  (6)  are  generic  control 
systems  judged  to  represent  useful  control  system  variations  for  an  experimen¬ 
tal  investigation  based  on  the  attack  helicopter  mission. 

In  general,  a  digital  representation  of  the  control  system  transfer  func¬ 
tions  is  obtained  by  the  use  of  the  E-transform;  using  programs  described  in 
reference  4,  the  appropriate  difference  equations  are  obtained  from  the 
corresponding  s-plane  transfer  functions. 

Block  diagrams  of  the  various  control  system  configurations  are  presented 
in  figure  i. 


«*Sfa 
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SCAS/ATTITUDE  HOLD: 


HOVER  AUGMENTATION  SYSTEMS: 


-1.75 


(b)  Roll. 

Figure  1.-  Continued. 
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MECHANICAL: 


LINKAGE 


DELC,  in. 

1  dfi  de9 

+ 

* 

in. 

1.0 


4 


THET0,  deg 


VERTICAL  AUGMENTATION  SYSTEMS: 


(d)  Collective. 
Figure  1.-  Concluded. 


To  provide  the  proper  control  force  characteristics,  the  simulator  cock¬ 
pit  is  equipped  with  a  control  loader  system  described  in  reference  4.  The 
reference  3  experiment  used  the  following  control  force  and  displacement  char¬ 
acteristics,  derived  from  reference  7,  to  simulate  the  advanced  helicopter. 


Controller 


Lateral 

Longitudinal 

Directional 

Collective 

Force  gradient,  lb/in. 

0.3 

0.5 

7.0 

Adjustable 

friction 

Breakout  force,  lb 

1.3 

1.0 

7.0 

Adjustable 

friction 

Travel,  in. 

±4.5 

±5.0 

±2.75 

0-12 

Mechanical  Flight  Controls 

The  baseline  mechanical  flight  control  system  uses  pilot  inputs  of 
(1)  longitudinal  cyclic  6e,  (2)  lateral  cyclic  6a,  (3)  directional  controls 
Sr,  and  (4)  collective  6C  to  determine,  respectively,  (1)  longitudinal 
swash-plate  angle  Bis,  (2)  lateral  swash-plate  angle  Ais,  (3)  tail  rotor 
collective  pitch  6TR,  and  (4)  main  rotor  collective  pitch,  0O.  The  relation¬ 
ships  between  the  pilot  control  position  and  control  surface  position,  derived 
from  reference  5,  are  as  follows. 

Longitudinal- 


Bis  =  5.0  -  3.0  <5e 


Lateral- 


Limits 


±5  in. 

+20°,  -10° 


Ai  =  -1.75  +  2.06  6„  Limits 
s  a 


6a: 

Als: 


Direct ional- 


±4.5  in. 
+7.5°,  -IV 


0^,R  =  9.25  -  8.45  &r  Limits  { 


Collective- 


90  =  1.0  +  1.46  <5C  Limits 


<5r:  ±2.75  in. 


0TR:  +32.5°,  -14' 


5C:  0-12  in. 


0O:  1-18.5' 
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Stability  and  Control  Augmentation  System  (SCAS) 


Limited  authority  SCAS  actuators  produce  additional  control  surface  motion 
in  response  to  sensed  aircraft  motion  parameters  (SAS)  and  pilot  control  inputs 
(CAS)  in  the  longitudinal,  lateral,  and  directional  axes.  The  SCAS  control 
mode  may  be  selected  by  the  researcher  for  each  axis  individually  or  for  all 
three  axes  collectively.  The  transfer  functions  for  the  SCAS  are  presented 
below  together  with  the  simplifications  employed  for  the  purposes  of  the 
simulation. 

Longitudinal  SCAS- 


where 


<5Bi 

«Bi  =  — r~ 


6B1  dBj. 

8  +  -  •  u  +  -z —  •  6e 

u  oe 


6B-l 

— 


8.54  s2(s  +  1.756)  10.62(s  +  0,3)(s  +  0.975) 

(s  +  0. 1)  (s  +  0.145)  (s  +  0.15) 


-  deg/rad 


Simplifying, 

.  a  8.54  s2(s  +  1.756)  +  10.62(s  +  0.1)(s  +  0.3)(s  +  0.975) 
6  ^s;  =  (s  +  0. 1)  (s  +  0.15) 

-  19-16(s3  +  1.545  s2  +  0.2327  s  +  0.01621) 

(s  +  0. 1) (s  +  0.15) 

19. 16(s  +  1 .386) [s2  +  2(0.72) (0. ll)s  +  (0.11)2] 

(s  +  0 . 1)  (s  +  0. 15) 


(s)  s  15. 71  (s  +  1.386) 


and 


6B 


l 


u 


(s) 


4.452  x  1Q~3  (s  +  0. 3)  (s  +  0.975) 
(s  +  0. 15) (s  +  1.0) 


deg/ft/sec 


Simplifying, 

-  (s)  =  8.681  x  io~3 

u 

Finally, 

661  .  -12.32  8(s  +  1.756) _ 

6e  ^s;  =  (s  +  0. 145) (s  +  0.147)(s  +  3.45) 
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Lateral  SCAS- 


6AX 


6A1  = 


6AX 

+  -z —  ■  6a 
6a 


where 


^  (s)  =  -1..461  s2(s  +.jL3)  _  1.45_s(s  +,2.28)  .  d e  /rad 
<P  (  )  (s  +  0.1)  (s  +  0.2)  (s  +  0.87)  aeg/raa 

=  -2.911  s (s  +  2. 3) (s  +  0.0175) (s  +  0.5686) 

(s  +  0 . 1) (s  +  0. 2) (s  +  0.87) 


Simplifying, 

6A 


1  (s)  a-  _i  90  s(s  +  2.3) 
(s)  -  i.w  (s  +  0  2) 


Finally, 


5Al  .  .  0.908  s (s  +  2.3) _ 

5a  (s  +  0.2)(s  +  0.2) (s  +  0.769) 


deg/in. 


Directional  SCAS- 


.  ^TR  >  ^TR 
6TR  r  *r  +  6, 


6r  + 


TR 


$  + 


TR 


where 


and 


5 

(s)  =  87,04  +Sq  ~  deg/rad/sec 


TR 


(s)  = 


-112  s 


(s  +  0. 2) (s  +  5.0) 


deg/in. 


For  V  <  50  knots,  6jR/4>  =  6TR/v  =  0.  For  V  >  5G  knots  only. 


TR 


TR 


(s)  = 


-324.3  Ks‘ 


614.8  Ks 


(s  +  0.  2)  (s  +  10)  (s  +  0. 2)  (s  +  10) 
where  K  *  0.5  -  0.00333(V  -  50)  (V  ~  knots).  Simplifying, 

_^TR  ,  .  _  -32.43  Ks (s  +  1.896) 

4,  ^s'  -  (s  +  0.2) 


~  deg/rad 
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and 


TR  ,  „  -831.4  ,  , 

~  (S)  =  V(7  +  14.7)  ~  ^g/ft/sec 


TR  ,  .  57.3 

—  (S)  S  "  ~ 


SCAS  limits-  SCAS  actuator  authority  limits  are  presented  in  reference  7 
as  percentages  of  equivalent  full  controller  deflection  as  follows: 

1.  10%  for  pitch  and  roll  SCAS 

2.  20%  for  yaw  SCAS 

The  following  control  surface  limits  result: 

6Bi  ->  ±3.0° 

<5Ai  -*  ±1.86° 

6tr  *  l9-3' 


Attitude  Hold 

According  to  reference  7,  the  attitude  hold  mode  is  available  below 
V  =  50  knots  by  switching  out  the  CAS  in  the  pitch  and  roll  axes,  that  is 

6Bj.  <5Bi 

6Bi  =  — —  •  6  +  -  •  u 

9  u 

and 

6Aa 

5Aj.  =  — —  •  <(» 

9 

and  by  providing  a  pseudo-heading-hold  feature  in  yaw,  that  is 

str  ‘  87-M  (i-roXrrM)  r 


Hover  Augmentation 

The  implementation  of  a  hover  position  hold  system  through  the  pitch  and 
roll  SCAS  actuators  consists  of  the  following  logic: 
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and 


6Al  =  K6a6a  +  K/fi  J  «a  +  KpP  +  +  K.yh  +  Kye^ 

where  the  h  subscript  indicates  positions  and  inertial  velocities  in  an  air¬ 
craft  heading-referenced  axis  system  with  origin  at  the  nominal  center  of 
gravity,  and  the  e  terms  indicate  position  errors  from  the  pilot-designated 
hover  point. 

The  standard  simulation  software  described  in  reference  4  calculates  the 
north  and  east  components  of  the  aircraft  inertial  velocity  (VN  and  VE, 
respectively).  The  transformation  from  these  earth-referenced  velocity  compo¬ 
nents  to  the  heading-referenced  components  (XDH  and  YDH)  utilizes  the  sine  and 
cosine  of  the  heading  angle  (SPSI  and  CPSI)  as  follows: 

XDH  =  VN*CPSI  +  VE*SPSI 

YDH  =  -VN*SPSI  +  VE*CPSI 

The  heading-referenced  position  errors  ex^  and  are  calculated 

through  an  integration  of  the  appropriate  velocity  components  which  commences 
when  the  pilot  designates  a  hover  point. 

These  heading-referenced  quantities  are  also  used  by  the  display  dynamics 
program  to  calculate  the  positions  of  various  symbols  on  the  pilot’s  electronic 
display. 


Heading  Hold 

With  the  heading  hold  mode  selected,  the  directional  axis  SCAS  equation 
becomes : 

STR  ‘  \Sr  +  K/5  J«r  +  M  +  V* 

The  intent  of  this  control  mode  is  to  provide  a  yaw  rate  command-heading  hold 
control  system  through  the  pilot's  directional  controls. 


Vertical  Augmentation 

With  vertical  augmentation  selected,  a  simulated  collective  SCAS  is 
implemented,  consisting  of  the  following  logic: 

\  ■  Vc +  K/«c  +  ^ +  Vh 
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The  objective  of  this  SCAS  mode  is  to  provide  an  altitude  rate  command- 
altitude  hold  control  system  through  the  pilot's  collective  stick. 


Summary  of  Equations 

In  general,  the  various  control  systems  to  be  investigated  are  implemented 
as  perturbations  on  the  basic  mechanical  flight  control  system;  that  is: 

A1S  =  -1.75  +  2.06*DELA  +  DELA1 

with  AlS  limited  to  +7.5°  to  -11°,  DELA  limited  to  ±4.5  in.,  and  DELAl  limited 
to  ±1.86°. 

BIS  =  5.0  -  3 . 0*DELE  +  DELB1 

with  BIS  limited  to  +20°  to  -10°,  DELE  limited  to  ±5  in.,  and  DELB1  limited  to 
±3.0°. 

THETTR  =  9.25  -  8.45*DELR  +  DELTR 

with  THETTR  limited  to  +32.5°  to  -14°,  DELR  limited  to  ±2.75  in.,  and  DELTR 
limited  to  ±9.3°. 

THET0  -  1.0  -f  1.46*DELC  +  DELTH0 

with  THET0  limited  to  1.0°  to  18.5°,  and  DELC  limited  to  0.0  to  12.0  in. 

The  perturbation  quantities  DELAl,  DELB1,  DELTR,  and  DELTH0  are  calcu¬ 
lated  using  logic  determined  by  the  control  mode  selected  (tables  1  through  4). 
(The  SCAS  actuator  limits  specified  above  are  nominal  values  and  may  be  set  to 
any  other  values  by  the  researcher.) 

TABLE  1.-  LOGIC 

_ in _ 

Control  mode  DELB1  = 

Pitch  SCAS  OFF  - 

Pitch  SCAS  ON  [ 7  a  a-^1*DELE  +  21.77*THETR 

L(s  +  0.145)(s  +  0.147)(s  +  3.45)J 

+  15 . 71*QB  +  0. 00868 1*UB 

Attitude  hold  21.77*THETR  +  15.71*QB  +  0.008681*UB 

Hover  hold  UKDELE*DELE  +  ^  *(UKDELEI.*DDELE  +  UKX*XDH)  +  UKTHETH*THETR 

+  UKQH*QB  +  UKXD*XDH 

where  DDELE  is  the  perturbation  of  DELE  from  its  value  at  the 
time  of  engagement  passed  through  a  dead  zone  of  ±0.1  in. _ 
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TABLE  2.- 


LOGIC 


Control  mode 

DELA1  = 

Roll  SCAS  OFF 

— 

Roll  SCAS  ON 

[(.  +  0°2)(“  +(0.2)(;31  0.769>]*DELA  -  l'Ks  +  oiD*88 

Attitude  hold 

-‘•9*C  +  o!3)*88 

Hover  hold 

UKDELA*DELA  +  j  * (UKDELAI*DDELA  +  UKY*YDH)  +  UKPHIH*PHIR 

+  UKPH*PB  +  UKYD*YDH 

where  DDELA  is  the  perturbation  of  DELA  from  its  value  at  the 
time  of  engagement  passed  through  a  dead  zone  of  ±0. 1  in. _ 


TABLE  3.-  6.  LOGIC 

0TR 

Control  mode 

DELTR  = 

Yaw  SCAS  OFF 

YA"  SCAS  0N  [-(S-+  o“)i°  I  5.0)],PELR  +  87 •  04*(l^X2)*RB 

+  KH*  £(-21.61  +  0.06391*UB)*^^^|^*PB  -  *VbJ 


Attitude  hold 
Heading  hold 


where 


UBKTS 


[c  7 ofta  ‘  5.o)]*°elr  +  87-04‘(nfe)(HT7r)*RB 
UKDELR*DELR  +  -  *(UKDELRI *DDELR  +  UKPSI*RB)  +  UKRH*RB 


where  DDELR  is  the  perturbation  of  DELR  from  its  value  at  the 
time  of  engagement  passed  through  a  dead  zone  of  ±0.1  in. _ 
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TABLE  4.-6.  LOGIC 

6o 

Control  mode 

DELTH0  = 

Collective 

SCAS  OFF 

— 

Altitude  hold 

UKDELC*DELC  +  -  * (UKDELCI*DDELC  +  UKAH*ALTD)  +  UKHD*ALTD 
s 

where  DDELC  is  the  perturbation  of  DELC  from  its 
time  of  engagement  passed  through  a  dead  zone  of 

value  at  the 
±0.1  in. 

4.  DISPLAY  DYNAMICS 


General 

The  purpose  of  the  display  dynamics  portion  of  the  mathematical  model  is 
to  produce  the  signals  used  to  drive  the  moving  symbols  on  the  various  candi¬ 
date  electronic  display  formats.  These  signals  are  either  simply  elements  of 
the  aircraft  state  vector  or  the  result  of  certain  logic  applied  to  selected 
state  vector  elements  to  produce  the  desired  dynamic  characteristics.  The 
moving  symbols  are  organized  in  this  section  on  the  basis  of  the  type  of  infor¬ 
mation  they  convey;  that  is,  orientation,  situation,  command,  and  fire  control. 

An  additional  function  of  this  portion  of  the  program  is  to  alter  the  dis¬ 
play  logic  as  a  function  of  four  discrete  display  modes  —  cruise,  transition, 
hover,  and  bob-up  —  which  are  selected  manually  by  the  pilot.  Reference  9 
describes  the  operational  requirements  associated  with  each  display  mode  as: 

1.  Cruise  —  high-speed  level  flight  en  route  to  the  forward  edge  of  the 
battle  area 

2.  Transition  —  low-speed  nap-of-the-earth  maneuvers,  such  as  dash,  quick 
stop,  and  sideward  flight 

3.  Hover  —  stable  hover  with  minimum  drift 

4.  Bob-Up  —  unmask  and  remask  maneuvers  over  a  selected  horizontal  ground 
position 

In  addition  to  the  electronic  display  symbol  drive  logic,  the  display 
dynamics  program  also  provides  signals  for  the  following  cockpit  instruments: 

1.  Attitude-director  indicator  (ADI) 

2.  Horizontal  situation  indicator  (HSI) 

3.  Radar  altimeter 
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4.  Barometric  altimeter 

5.  Instantaneous  vertical  speed  indicator  (IVSI) 

6.  Normal  accelerometer 

7.  Airspeed  indicator 

8.  Engine  torque 


Flight  Control  Display  Logic 

A  typical  electronic  display  format  is  illustrated  in  figures  2  through  5. 
The  primary  symbols  used  by  the  pilot  to  control  his  aircraft  are  the  velocity 
vector,  cyclic  director  symbol,  and  hover  position  symbol  (fig.  3).  The  logic 
and  scaling  of  the  parameters  that  drive  these  symbols  vary  as  a  function  of 
display  mode. 

Cruise  mode-  These  symbols  do  not  appear  in  the  cruise  display  mode. 

Transition  mode-  The  velocity  vector  is  driven  directly  by  the  horizontal 
components  of  Doppler  velocity  in  the  transition  mode;  that  is,  displayed 
vertical  motions  of  the  vector  are  driven  by  the  longitudinal  component  of 
heading-referenced  velocity  (XDH)  while  its  lateral  motions  are  driven  by  YDH. 

The  displayed  vertical  motion  of  the  cyclic  director  symbol  with  respect 
to  the  tip  of  the  velocity  vector  is  driven  by  washed-out  pitch  attitude  with 
a  washout  time  constant  of  50  sec.  Laterally,  the  symbol  is  driven  by  roll 
attitude  for  roll  angles  greater  than  5.73°  and  by  washed-out  roll  attitude 
for  smaller  values  of  roll  angle.  For  the  latter  case,  the  washout  time  con¬ 
stant  is  10  sec. 

Hover  mode-  For  the  smaller  values  of  velocity  encountered  in  the  hover, 
the  velocity  vector  is  driven  by  a  complementary-f iltered  estimate  of  the 
longitudinal  and  lateral  heading-referenced  inertial  velocities,  x  and  y, 
respectively. 

These  estimates  are  obtained  by  blending  low-frequency  Doppler  velocity 
data  with  high-frequency  inertial  acceleration  data.  Assuming  1-g  flight  and 
small  values  of  pitch  and  roll  attitudes,  the  inertial  accelerations  may  be 
expressed  as: 


x  s  aXcg  “  8e 
y  s  ayCg  + 


The  velocity  estimates  are  expressed  as: 


(txs  +  l)X  +  [(s  +  Ax)(s  +  1/Tx)](3xcg  "  g6) 
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Figure  2.-  Electronic  display  format  (ref.  8). 
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AIRCRAFT  REFERENCE 
SYMBOL 


VELOCITY  VECTOR 


O - CYCLIC  DIRECTOR 

SYMBOL 


HOVER  POSITION 
SYMBOL 


HORIZON  LINE 


Figure  3.-  Central  symbology. 
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Figure  4.-  Peripheral  symbology. 
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Figure  5.-  Fire  control  symbology. 
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(v  +  l)y  +  {(s  + Ay)[s+  (1/Ty)]}(aycg  +  g^) 


The  filter  time  constants  are  retained  as  variables, 
suggested  in  reference  9  is  to  assume  that  both  ax 


A  further  simplification 

and  av  are  negligible 
eg  yeg  & 

with  respect  to  the  gravitational  acceleration  components.  Both  the  complete 

and  the  simplified  expressions  for  x  and  y  are  implemented  in  the  program. 


The  cyclic  director  symbol  is  driven  by  washed-out  pitch  attitude 
(10-sec  time  constant)  and  washed-out  roll  attitude  (10-sec  time  constant). 


These  changes  in  logic  occur  instantaneously  at  the  time  of  the  switch 
from  transition  to  hover  mode. 

Bob-up  mode-  The  logic  driving  the  velocity  vector  and  cyclic  director 
symbol  remains  the  same  as  the  hover  mode  logic.  The  hover  position  symbol  is 
now  driven  vertically  by  EXH  and  laterally  by  EYH  where  EXH  and  EYH  are  the 
integrals  of  XDH  and  YDH,  respectively,  with  integration  commencing  at  the 
time  the  bob-up  display  mode  is  selected.  Finally,  a  command  heading  symbol, 
which  has  remained  fixed  on  the  display,  is  now  driven  by  the  difference 
between  the  current  heading  and  the  heading  that  existed  at  the  time  the  bob-up 
display  mode  was  selected. 


Torque  Equation 

One  of  the  status  parameters  displayed  to  the  pilot  is  engine  torque  in 
percent  (fig.  4).  Since  no  engine  or  rotor  dynamics  are  considered  in  the 
model,  an  approximate  relationship  between  engine  torque  and  main  rotor  col¬ 
lective  pitch  is  used  to  drive  this  symbol.  This  relationship  is  derived  as 
follows.  The  relationship  between  power  required  and  power  turbine  torque  may 
be  expressed  as: 


where 


P  =  shaft  power  =  power  required,  (ft-lb)/sec 
N2  =  power  turbine  speed  =  2*104  rpm 
Q  =  power  turbine  torque,  ft-lb 

In  hover,  the  total  power  coefficient  is  equal  to  the  sum  of  the  induced,  pro¬ 
file,  and  miscellaneous  power  coefficients: 


"Pi 


+  C„  +  C 


Pmisc 
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Typically 


1.67  CL 


where 


For  the  advanced  helicopter  in 


P  = 


hover  (CT  =  0.0066) 
6.35  x  10~4 

1.05  x  ios  (f t-lb) /sec 


HP 


P 

550 


1900 


This  value  of  hover  horsepower  required  is  approximately  62%  of  the 
available  3  086  hp  in  the  advanced  helicopter.  Therefore  the  torque  required 
to  hover  out  of  ground  effect  is  assumed  to  be  62%.  To  determine  the  relation¬ 
ship  between  a  perturbation  from  the  trim  value  of  main  rotor  collective  pitch 
and  the  resultant  perturbation  in  torque  from  the  nominal  hover  value,  the 
following  technique  is  used.  Since  cp  =  cPi  +  cPo  +  Cpmisc, 


3Cp  3CPi  9Ct 
39q  3C-j  30q 


and 


Therefore, 


but 


and 


OR 


3Q  _  _3P_  60 

30o  36q  ‘  2 it 


3Cpi  #  _3T_ 
90q 


RR  •  •—  N 


2 
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Therefore, 


For  the  advanced  helicopter  in  hover  out  of  ground  effect: 

Cip  =  0.0066 

m  =  453.42  slugs 

Zo  =  -4.562  ft/sec2/deg 
uo 

OR  =  (30.26  rad/sec) (24  ft) 

N2  =  2  x  104  rpm 

Therefore, 

^  4  61.86  (f t-lb) /deg 

The  maximum  torque  in  foot-pounds  is  estimated  by  dividing  the  hover  value 
of  torque  in  foot-pounds  by  that  value  in  percent;  that  is, 

Q  -  49?:]s  =  811.6  ft-lb 
max  0.6157 


Therefore , 


3Q  i  61.86 
360  811.6 


100  =  7.6%/deg 


This  hover  approximation  for  torque  is  also  used  for  forward  flight  since  the 
experiment  for  which  the  model  was  designed  (ref.  3)  only  involved  hover  and 
low-speed  flight. 


Fire  Control  Display  Logic 

Figure  5  illustrates  typical  symbology  used  to  assist  the  pilot  in  his 
weapon  delivery  task.  The  logic  driving  the  cued  line  of  sight  symbol  causes 
the  symbol  to  overlay  a  simulated  fixed  ground  target  in  the  background  image. 
The  background  imagery  simulates  the  video  received  from  a  sensor,  such  as 
forward-looking  infrared  or  low-light-level  television  mounted  at  some  loca¬ 
tion  on  the  aircraft.  This  location  is  variable  to  allow  the  investigation  of 
the  effects  caused  by  the  sensor  offset  from  the  pilot's  eye  position. 

The  derivation  of  the  logic  for  the  cued  line  of  sight  symbol  proceeds 
as  follows.  Let  (xT,  y^,  z^)  represent  the  target  position  in  an  aircraft  body 
axis  system  with  origin  at  the  sensor  location.  Therefore  the  desired  values 
of  displayed  target  azimuth  and  elevation  are: 
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-1  yT 
Az  =  tan  — 
XT 


E«  -  -sl"_1  rs 


where  Rs  =  /xT2  +  yT2  +  z^2 .  The  values  of  xT,  yT,  and  zT  are  calculated  as 
follows.  Target  position  in  an  earth-referenced  coordinate  system  with  origin 
at  the  sensor  location  is  expressed  as: 

xTp  =  R±  cos(Az.  +  «,.)  -  Ax 
yTp  «  R±  sin(Azi  +  <p±)  -  Ay 

ZTP  =  Zs  hTER 


where 


R^  =  initial  target  range,  ft 

A^  =  initial  target  azimuth,  deg;  positive  right 

4^  ■  aircraft  heading  at  time  of  target  designation,  deg 

Ax  =  translational  motion  of  sensor  in  N-S  direction  after  target  designa¬ 
tion,  ft;  positive  north 

Ay  =  translational  motion  of  sensor  in  E-W  direction  after  target  designa¬ 
tion,  ft;  positive  east 


zs  =  simulated  sensor  height  above  terrain  board  reference,  ft 

h_„D  =  mean  simulated  height  of  terrain  above  terrain  board  reference,  ft 
TEK 

Therefore,  transforming  to  a  body  axis  system  with  origin  at  the  sensor 
location: 


cos  6  cos  ip 


cos  9  sin  \p 


-sin  0 


sin  <j>  sin  9  cos  \p  sin  <j>  sin  0  sin  ip  sin  <p  cos  0 

-  cos  <ji  sin  \p  +  cos  <p  cos  ip 

cos  <j>  sin  0  cos  ip  cos  sin  ®  sin  \p  cos  <{>  cos  0 

|_  +  sin  <f>  sin  ip  -  sin  <f>  cos  ip 


XTP 

yTP 

ZTP 
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Summary  of  Equations 


... - 


Orientation-  The  following  parameters  are  used  to  drive  the  moving  symbols 
which  provide  information  on  aircraft  orientation: 

Symbol  Parameter 

AIRCRAFT  HEADING  PSI 

HORIZON  LINE  THET,  PHI 


Situation-  Aircraft  position  and  velocity  information  in  the  horizontal 
and  vertical  planes  are  provided  to  the  pilot  through  the  following  symbols: 


Symbol 


Parameter 


VELOCITY  VECTOR 


HORIZONTAL 


(LONGITUDINAL  AIRSPEED 
HOVER  POSITION 


XDH,  YDH  (TRANSITION) 

XDHAT,  YD HAT  (HOVER/BOB-UP) 

UBKTS 

EXH,  EYH  (BOB-UP) 


VERTICAL 


RADAR  ALTITUDE 
RATE  OF  CLIMB 


HAGL 

ALTD 


The  velocity  vector  symbol  is  driven  in  the  transition  mode  by  the  true 
values  of  ground  velocity,  XDH  and  YDH,  and  in  the  hover  and  bob-up  modes  by 
estimates  of  these  values,  XDHAT  and  YDHAT,  respectively.  These  estimates  are 
derived  by  complementary  filtering  techniques  described  in  section  4  (Flight 
Control  Display  Logic) .  The  ability  to  vary  the  scaling  of  the  velocity  vector 
is  retained  in  the  display  dynamics  program.  Thus: 


IUKDXD*XDH  (TRANSITION) 
UKDXD*XDHAT  (HOVER/BOB- UP) 

and 

(UKDYD*YDH  (TRANSITION) 
UKDYD*YDHAT  (HOVER/BOB -UP) 


where  UKDXD  and  UKDYD  are  constants  the  values  of  which  may  be  selected  by  the 
researcher  and  which,  in  general,  vary  as  a  function  of  display  mode. 

In  the  bob-up  mode,  the  hover  position  symbol  moves  in  response  to  the 
variables  EXH  and  EYH.  Thus: 


HQVX  =  UKDX*EXH 
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and 


HOVY  =  UKDY*EYH 

where  UKDX  and  UKDY  are  constants  whose  values  may  be  selected  by  the 
researcher. 

Additional  status  information  includes  engine  torque  and  lateral  acceler¬ 
ation.  The  expression  for  engine  torque  is  derived  in  section  4  (Torque 
Equation).  The  torque  response  to  collective  pitch  is  lagged  by  a  first-order 
filter  with  a  0.1-sec  time  constant.  Thus: 


TRQ  =  TRQTRM  +  (THET0  -  TRQCB) *TRQG* 


10 

s  +  10 


where 


TRQTRM  =  trim  value  of  percent  torque  in  hover 
TRQCB  =  trim  value  of  collective  pitch  in  hover 
TRQG  =  percent  torque  change  per  degree  of  collective  pitch 


Lateral  acceleration  is  driven  by  the  parameter  AYP. 

Command-  The  cyclic  director  symbol  provides  "command"  information  in  the 
horizontal  plane  which,  if  properly  designed,  allows  the  pilot  to  reach  and 
maintain  a  stable  hover.  Thus, 


VTIPX  =  WECX  +  UKDTHT*THET* 


VTIPY  =  WECY  +  UKDPHI*PHI 


where  UKDTHT  and  UKDPHI  are  constants  the  values  of  which  may  be  selected  by 
the  researcher  and  which,  in  general,  vary  as  a  function  of  display  mode;  the 
nominal  values  of  Ti  and  T2  are  functions  of  display  mode  as  follows: 

Transition  Hover /bob-up 

Ti,  sec  50  10 

T2 ,  sec  10  for  PHIR  <0.1  10 

«=  for  PHIR  >  0. 1 

In  addition,  a  command  heading  symbol  is  provided;  this  symbol  is  driven  by 
the  difference  between  the  current  heading  and  the  heading  that  existed  at  the 
time  the  bob-up  display  mode  was  selected  (EPSIBU). 
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Finally,  logic  for  a  collective  stick  director  is  provided.  The  director 
logic  is  implemented  as  a  weighted  sum  of  altitude  and  altitude  rate  which 
drives  the  original  rate  of  climb  symbol;  thus, 

ALTDRC  +  UKDALTD  *ALTD  +  UKDHAGL* (HAGL-100) 

For  rate  of  climb  information  only,  UKDHAGL  is  set  to  zero. 

Fire  control-  The  equations  derived  in  section  4  (Fire  Control  Display 
Logic)  for  target  azimuth  and  elevation  angles,  TADSX  and  TADSY,  respectively, 
are  implemented  as 

XTP  =  RGT*COS (AZAS*D2R)  +  DELX 
YTP  =  RGT*SIN(AZAS*D2R)  -  DELY 
ZTP  =  HPR  -  HTER 

where 

AZAS  =  AZA  +  PSIGT 
DELX  =  YPR  -  YPR0 
DELY  =  XPR  -  XPR0 

and  RGT  and  AZA  are  constants  selected  by  the  researcher  to  represent  initial 
target  range  and  azimuth,  respectively. 

XT  =  Til *XTP  +  T12*YTP  +  T13*ZTP 

YT  =  T21*XTP  +  T22*YTP  +  T23*ZTP 

ZT  =  T31*XTP  +  T32*YTP  +  T33*ZTP 

TADSX  =  R2D*ATAN2(YT,XT) 

TADSY  =  -R2D*ASIN(ZT/SLANTR) 
where  SLANTR  =  /XT2  +  YT2  +  ZT2. 

5.  WIND  AND  TURBULENCE 

General 

As  described  in  reference  4,  a  model  of  atmospheric  turbulence  (BWIND) 
based  on  MIL-F-8785B  (ref.  10)  is  implemented  in  the  standard  simulation  soft¬ 
ware.  Inputs  required  for  the  model  are  airspeed  (VRW) ,  the  scale  lengths 
(UAL,  VAL,  WAL),  and  the  turbulence  intensities  (UDISP,  VDISP,  WDISP) .  Outputs 


of  the  standard  program  are  linear  and  angular  turbulence  velocities  in  body 
axes  (UTURB,  VTURB,  WTURB,  PTURB,  QTURB,  RTURB).  A  suggested  revision  to 
MIL-F-8785B  (ref.  11)  contains  a  low-altitude  turbulence  model;  the  changes 
basically  involve: 

1.  Different  relationships  between  scale  lengths  and  altitude 

2.  Different  relationship  of  turbulence  intensities  to  mean  wind  speed 
and  altitude 

3.  Linear  turbulence  velocities  in  wind  axes  rather  than  body  axes 

Approximations  to  these  changes  are  implemented  in  the  turbulence  model  for 
the  advanced  helicopter. 

In  addition  a  model  of  the  steady  wind,  which  includes  provisions  for 
wind  shear  in  magnitude  and  direction,  is  included  in  the  program. 

Steady  Wind 

To  provide  the  effects  of  steady  wind  and  wind  shear,  the  magnitude  of 
the  steady  wind  (Vw)  is  specified  at  two  altitudes:  20  ft  and  200  ft,  AGL. 
Linear  interpolation  is  used  to  determine  mean  wind  speed  between  these  alti¬ 
tudes.  Beyond  these  altitude  extremes  the  mean  wind  speed  remains  constant. 
Wind  direction  may  be  specified  as  a  function  of  altitude  in  a  similar  fashion. 
The  components  of  the  steady  wind  in  an  earth-axis  system  are  then  derived  for 
use  as  inputs  to  the  standard  simulation  software. 


Turbulence 


The  scale  lengths  required 
presented  in  reference  11,  that 


as  inputs  to  BWIND  are  approximations  to  those 
is, 


and 


^w  = 


Lu  =  Lv  = 


for 

{20 

for 

f 

5h 

for 

100 

for 

1000 

for 

h  >  20  ft 
h  <  20  ft 

200  ft  >  h  >  20  ft 
h  <  20  ft 
h  >  200  ft 


To  provide  a  valid  turbulence  model  when  airspeed  is  zero,  steady  wind 
speed  (Vw)  is  used  as  the  minimum  value  for  the  required  airspeed  input. 

Finally,  the  vertical  turbulence  intensity  is  specified  in  reference  11 
as  being  10%  of  the  mean  wind  speed  measured  at  20  ft,  AGL.  The  ratio  of  the 
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horizontal  turbulence  intensities  to  the  vertical  intensity  varies  as  a  func¬ 
tion  of  altitude  from  a  value  of  1.0  at  1,000  ft  to  2.0  at  zero  altitude.  A 
value  of  2.0  is  selected  as  being  appropriate  for  the  purposes  of  any  simula¬ 
tions  conducted  using  this  model. 

The  outputs  of  BWIND  are  normally  the  six  turbulence  velocities  in  body 
axes.  Two  revisions  to  the  output  velocities  are  made  for  this  model: 

1.  The  angular  turbulence  velocity  components  are  set  to  zero. 

2.  The  linear  turbulence  velocity  components  are  assumed  to  be  in  wind 
axes  and  are  then  transformed  to  earth  axes  and  summed  with  the  steady  wind 
components . 

In  the  frequency  domain,  the  high-frequency  cutoff  of  the  turbulence  com¬ 
ponents  occurs  in  general  at  a  frequency  of  V/L,  according  to  reference  10. 

As  an  example,  assume  the  mean  wind  speed  to  be  constant  at  15  knots  and  allow 
the  altitude  to  vary  from  20  to  200  ft  AGL.  The  following  variations  in 
turbulence  break  frequencies  result  for  hovering  flight: 

Altitude  Vw/Lw,  rad /sec  Vw/Lu  =  Vw/Lv,  rad /sec 


20  ft  1.27  0.25 

200  ft  0.13  0.025 


Summary  of  Equations 

Steady  wind-  The  magnitude  of  the  steady  wind  is  specified  at  two  values 
of  HAGL;  linear  interpolation  is  used  to  determine  wind  speed  between  these 
values  of  altitude,  that  is, 


Wind  direction  (PSIW)  is  then  specified  as  a  constant  or  a  function  of 
altitude.  The  components  of  the  steady  wind  in  earth  axes  are  then  calculated 


VNW  «  -VW*COS(PSIW) 


VEW  =  -VW*S1N(PSIW) 

VDW  =  0 

Turbulence-  The  scale  lengths  required  as  inputs  by  BWIND  are  specified  as 

(HAGL  for  HAGL  >  20 
WAL  =  < 


HAGL 

for 

HAGL  >  20 

20 

for 

HAGL  <  20 

5*HAGL 

for 

200  >  HAGL 

> 

20 

100 

for 

HAGL 

< 

20 

1000 

for 

HAGL 

> 

200 

The  minimum  value  of  the  BWIND  input  velocity  VRW  is  set  equal  to  the 
steady  wind  speed  VW  and  the  turbulence  intensities  are  specified  as  follows: 

WDISP  =  0.1*VWI 

UDISP  =  VDISP  =  0. 2*VWI 

The  three  turbulence  velocity  components  that  are  outputs  of  the  program 
are  transformed  from  wind  axes  to  earth  axes  as: 

VNTURB  =  -UTURB*C0S (PSIW)  +  VTURB*SIN(PSIW) 

VETURB  =  -UTURB*SIN (PSIW)  -  VTURB*COS (PSIW) 

VDTURB  =  WTURB 

Finally  the  earth-axis  components  of  the  steady  wind  are  summed  with  these 
turbulence  velocity  components,  yielding, 

VTWN  =  VNW  +  VNTURB 

VTWE  =  VEW  +  VETURB 

VDW  =  VDTURB 

6 .  CONCLUDING  REMARKS 


This  report  has  described  the  derivation  of  a  simplified  mathematical 
model  of  a  helicopter  suitable  for  use  in  piloted  simulator  experiments  involv¬ 
ing  tasks  that  may  include  relatively  large  variations  in  airspeed.  The  model 
was  designed  for  a  particular  program  of  research  which  investigated  the 
effects  of  control  system  and  display  variations  for  an  attack  helicopter  mis¬ 
sion;  however,  with  appropriate  modifications  to  the  trim,  stability,  and 
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control  parameters,  it  may  also  be  used  for  the  simulation  of  other  types  of 
helicopters  performing  different  tasks.  The  major  limitations  of  the  model 
are:  (1)  only  six  degrees  of  freedom  are  simulated  —  no  engine  or  rotor 

dynamics  is  included;  and  (2)  for  a  given  airspeed,  the  aerodynamic  forces  and 
moments  are  expressed  as  linear  approximations  of  the  actual  forces  and 
moments. 


APPENDIX 


AERODYNAMIC  PARAMETERS 


The  appendix  contains  a  summary  of  the  parameters  which  characterize  the 
aerodynamic  forces  and  moments  of  the  simulated  helicopter.  Table  A-l  contains 
the  mass  and  geometry  constants  required  as  inputs  to  the  program.  Table  A-2 
presents  the  reference  values  of  the  aircraft  state  and  control  parameters  for 
level  flight  as  a  function  of  airspeed.  Table  A-3  contains  the  nonzero  values 
of  the  aerodynamic  forces  and  moments  along  the  reference  trajectory. 

Tables  A-4  through  A-9  present  the  values  of  the  model's  stability  and  control 
parameters  as  a  function  of  airspeed. 


TABLE  Ai.-  MASS  AND  GEOMETRY  CONSTANTS 


Programming 

symbol 

Engineering 

symbol 

Definition 

Units 

Nominal 

value 

XI XX 

I 

5,140 

XX 

Body  axis  moments  of 

XIYY 

Iyy 

■ 

inertia 

slug-f t 

31,850 

XIZZ 

^zz 

30,850 

XIXZ 

*xz 

Cross-product  of  inertia 

slug-ft2 

1,260 

XMASS 

ra 

Aircraft  mass 

slugs 

453.42 

XP 

f  Pilot's  design  eye 

5.12 

YP 

Position  in  body 

ft 

.0 

ZP 

[  Axis  coordinates 

-2.39 

XPNVS 

Sensor  location  in 

15.5 

YPN  VS 

body  axis 

ft 

.82 

ZPNVS 

coordinates 

.0 
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TABLE  A2.-  REFERENCE  TRAJECTORY  PARAMETERS 


UBKTS 

Programming 

symbol 

Engineering 

symbol 

-40 

0 

40 

80 

120 

160 

Units 

WR 

WR 

0.0 

0.0 

6.66 

10.8 

4.46 

-1.89 

ft/sec 

A1SR 

Als 

SR 

-.16 

-.16 

-.44 

-.61 

-.90 

-1.19 

deg 

B1SR 

Bls 

SR 

-.45 

-.45 

2.97 

4.86 

7.17 

9.47 

deg 

THET0R 

CD 

0 

15.75 

15.75 

14.03 

13.23 

15.09 

16.95 

deg 

THETTRR 

9TRr 

21.46 

21.46 

13.44 

8.69 

8.02 

7.35 

deg 
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TABLE  A3.-  REFERENCE  VALUES  OF  FORCES  AND  MOMENTS 


UBKTS 

Programming 

symbol 

Engineering 

symbol 

-40 

-20 

0 

20 

40 

60 

80 

100 

120 

140 

160 

Units 

XKEF 

3.539 

2.934 

2.308 

2.59 

2.806 

2.692 

2.151 

1.503 

0.4095 

-0.684 

-2.222 

ft /sec' 

YREF 

yr 

2.047 

2.049 

2.050 

1.710 

1.365 

1.187 

1.009 

1.203 

1.396 

1.587 

1.777 

ft /sec' 

ZREF 

ZK 

-30.74 

-31.39 

-32.05 

-30.68 

-30.68 

-31.98 

-31.27 

-32.38 

-33.08 

-32.67 

-31.49 

ft /sec" 

UMREF 

"r 

-.169 

-.0845 

.0 

-.0473 

-.0947 

-.190 

-.2842 

-.2845 

-.2685 

-.2525 

-.2203 

r.id  /  sec ' 
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TABLE  A4 . -  X-FORCE  STABILITY  AND  CONTROL  PARAMETERS 
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TABLE  A5.-  Y- FORCE  STABILITY  AND  CONTROL  PARAMETERS 


TABLE  A6.-  Z-FORCE  STABILITY  AND  CONTROL  PARAMETERS 
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TABLE  A8.-  M-MOMENT  STABILITY  AND  CONTROL  PARAMETERS 
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TABLE  A9.-  N-MOMENT  STABILITY  AND  CONTROL  PARAMETERS 
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